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Abstract 

The relation between the distribution of work performed on a classical system by an external 
force switched on an arbitrary timescale, and the corresponding equilibrium free energy difference, 
is generalized to quantum systems. Using the adiabatic representation, we show that this relation 
holds for isolated systems as well as for systems coupled to a bath described by a master equation. 
A close formal analogy is established between the present "classical trajectory" picture over popu- 
lations of adiabatic states and phase fluctuations (dephasing) of a quantum coherence in spectral 
lineshapes, described by the stochastic Liouville equation. 
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I. INTRODUCTION 

Jarzynski had established a remarkably simple and general relationship between 

the distribution of work performed on a classical system by an external force and the free 
energy difference between the initial and final states . This interesting prediction, 

connecting equilibrium quantities (free energies) to nonequilibrium trajectories correspond- 
ing to measurements performed on an arbitrary timescale, has recently been verified in single 
molecule measurements M. 

n 

In this letter we utilize the adiabatic representation [7( to prove that this relation holds 
equally for quantum systems. For an isolated driven quantum system (i.e., not coupled to 
a bath) the work does not depend on the path, and the proof follows directly from a sum 
rule of nonadiabatic couplings. We further develop a path integral representation of systems 
coupled to a bath described by a master equation in the adiabatic basis using "trajectories" 
over adiabatic state populations, and generalize Jarzynski's relation to obtain the joint 
distribution of work and energy in a nonequillibrium measurement. A close resemblance 
is established between the physical picture underlying the distribution of work and the 
distribution of phase acquired by a quantum coherence due to coupling to a bath. 

Consider a quantum system coupled to an external force which changes its Hamiltonian 
H{t) by switching a parameter X{t) from A = Oatt = OtoA = lat time t. Initially 
the system is in thermal equilibrium and the instantaneous energy levels (eigenvalues of 
H{t)) are denoted e„(t) with the corresponding eigenstates |y9„(t)). The canonical partition 
function for a given value of A is 

Zx = J2 [-P^n{t)] = exp{-(3Fx) (1) 

n 

where Fx is the Helmholtz free energy. 

Let us consider an externally driven, but otherwise isolated, system. In the adiabatic 
representation 0| we expand its density matrix as 

pit) = ^ Pnm{t)\^n{t) >< ^mit)\, (2) 
nm 

where the coefficients Pnm{t) satisfy the Liouville equation 

d 



dt 



Pkl{t) = -iuJkl{t)pkl{t) - ^ Skl,mn{t)Pmn{t) (3) 



with ujkiit) = ekit) - eiit) and 



Skl,mn{t) = {'^k{t)\'^m{t))6in + {'^n{t)\(pi{t))6km 



(4) 



The nonadiabatic couphng S is of purely quantum origin and has no classical analogue. 

For a given realization of A(t), the solution of Eq. Q allows us to compute the conditional 
probability Kmn{t) of the system to be in the state \(prn{t)) at time t given that it started at 
state \ipn{0)), so that Pmm(t) = Kmn(t)pnn{0)- This probability is normalized by summing 
over final states J2m ^mnif) = 1. A second, less obvious, sum rule is obtained by summation 
over initial states K^nit) = 1- When the external force varies slowly, it follows from the 
adiabatic theorem [7[ that no nonadiabatic transitions take place, Kmn{t) = Smn, and this 
sum rule trivially holds. For finite switching timescales there are nonadiabatic transitions 
and the instantaneous eigenstates \ipn(t)) are no longer the solutions of the time dependent 
Schrodinger equation. Nevertheless, the second sum rule follows from the relation 



A direct consequence of Eq. (0) is that the following uniform distribution in adiabatic 
population space, obtained by setting Pnmit) = ^nm at all times in Eq. (j21) 



satisfies the Liouville equation. This interesting identity gives rise to the second sum rule, 
which holds despite the nonadiabatic transitions and is the key for proving the quantum 
Jarzynski identity. Alternatively, Eq. © can be viewed as a unitary transformation of the 
unit operator between adiabatic basis sets at different times. This equation then simply 
states that the unit operator is invariant to a unitary transformation. 

We shall now compute the ensemble average of exp {—jSW) where W is the cumulative 
work performed up to time t. If the system starts in state |<^„(0)) and ends up at time t in 
the state \ipm(t)), then the work made by the external force is W = emit) — e„(0), and we 
have 




(5) 



n 




(6) 



n 



{exp i-pW)) = —\exp [-Pen{0)]Kmn{t)exp [-/?(e„(t) - e„(0))]. 
^0 — 

ran 



(7) 



This gives 
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Using the second sum rule we can carry out the n summation and obtain Jarzynski relation, 



{exp{-pW)) = Z,/Zo (9) 

where Zi is the final state partition function for A = 1. Eq. (jH)) implies that it is generally 
possible to compute the free energy change Fi — Fq = —j3'^ \og{Zi/Zo) using the distribution 
of work. 

We next turn to a system coupled to a bath and described by a master equation in the 
adiabatic basis . We denote the transition rate from state n' to n by Rnn'if)- R depends 
on t since it is recast in the time dependent, adiabatic, basis. It further satisfies the detailed 
balance condition 

Rnn'(t)/Rn'n(t) = CXp {- (3uJnn' {t)) (10) 

For an isolated system the work equals the energy change and is therefore a state function. 
When the system is not isolated, the work becomes path dependent and may not be com- 
puted from the initial and final states alone, as was done in Eq. (jZj). For our model, the 
change in energy along the path stems from the variation of the eigenvalues e„(t) which is 
induced by the driving force, as well as the changes in their occupations, which are induced 
by the coupling to the bath (Fig. 1). The former is the work whereas the latter is the heat. 
Both work and heat are path dependent whereas the overall energy change em{t) — €^^(0) is 
a state function. 

We shall compute the joint distribution of work performed iW) and the change in system's 
energy {E) in the process by defining the joint energy-work generating function 

5(7, 5- 1) = {exp [-^W{t) - bE{t)\), (11) 

where (....) denotes an ensemble average over trajectories for a fixed realization of A(t). The 
moments of W and E can then be computed as derivatives with respect to the parameters 
7 and b 

{W{tYE\t)) = (-1)^^'^^^5(7, 5; t) U=s=o (12) 

By expanding the solution of the master equation perturbatively in the off diagonal 
elements of Rnn' {n 7^ n') it can be represented as a path integral in the adiabatic population 
space (see Fig. 1): Consider a family of trajectories which start at state |y9„(0)) at t = and 
end in state \ipjn{t)) at time t. In the course of a given trajectory the system had jumped 
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j times between different adiabatic states witli j = 0,1,2 Tlie conditional probability 

of the system to start at state |v9„(0)) and end in the state \(pm{t)) is then given by a sum 
over all possible values of j. The work can be computed by integrating den/dt along the 
continuous segments of each trajectory. The contribution of this family to the generating 
function can be calculated by solving the equation of motion for the following Green function 

= Yl Rmn'it)GUt) - l^GUt) (13) 
n' 

with the initial condition G'mn(O) = 6mn- It immediately follows from the path integral 
representation of the solution of the master equation that the generating function is given 
by averaging the Green function over the equilibrium distribution of initial states and a 
weighted summation over final states 

5(7, S-t) = y2 exp[-6[eUt) - en(0)]]G^„(t; ^^^^EtplM, (14) 

nm 

Eq. (fT^ provides the following physical picture for the process: We consider an ensemble 
of trajectories over adiabatic state populations. In the classical case |l| the ensemble is 
generated by the distribution of initial states, whereas here the ensemble originates from 
the random nature of the stochastic jumps among adiabatic states. The Green function 
represents an open system where population can flow to (or from) state m with a rate 
—'-fdem/dt which depends on the variation of the adiabatic energies with time and on the 
parameter 7. As a result, the total population '^^Gmnit) is not conserved and its value 
determines the generating function for the work and energy. The adiabatic representation 
thus provides a "classical trajectory" picture for the joint distribution of work and energy. 
To recover the Jarzynski relation we set 7 = /9 and 5 = in Eq. (jl4j) which then assumes 
the form 

S{(3;0,t) = J2Gm{f^,t) (15) 



where 



GMt) = J2GUP,tf-^^Ptp^ (16) 



Let us consider the (unnormalized) distribution of states m 



G„(ftt) = ££?<^f=iM (17) 

^0 
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It can be easily verified that Eq. (fTTj) is a solution to Eq. (fT^ : this follows from the detailed 
balance of R that guarantees that the first term in the r.h.s. of Eq. (fT^ vanishes when 
acting on Eq. (fT7|). Substituting Eq. (fTTj) in Eq. (fT^ immediately gives the Jarzynski 
relation (Eq. Q). It will be interesting to extend these results to systems described by 
more elaborate kinetic schemes such as fractional kinetics |,9| rather than the ordinary master 
equation. 

The validity of the Jarzynski relation for a system coupled to a bath at constant tem- 
perature has been estabhshed by Jarzynski Q and Crooks ^. The present derivation and 
the definition of work closely resemble Jarzynski's master equation approach. However, the 
path integral picture of continuous segments with discrete jumps provides a new insight 
for quantum systems and helps establish an interesting analogy with the theory of pure 
dephasing of quantum coherence, as will be shown next. 
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12 1, which has 



Eq. ()13|) closely resembles the stochastic Liouville equation of Kubo 
been widely used in the theory of NMR and optical lineshapes as well as in single molecule 
spectroscopy Computing the work generating function is then formally equivalent to 
computing the dephasing of a two level system {g and e) coupled to a bath. In that problem 
we consider equilibrium fluctuations of a quantum coherence a between the two levels which 
are coupled to a classical bath represented by a phase space F. a satisfies the stochastic 
Liouville equation 1^ 

= J2 WiT, T')aiT\ t) - ^UiT)aiT, t). (18) 
r' 

Here is a Markovian rate matrix and U (F) is the fluctuating energy gap between the two 
levels, a serves as a generating function for the absorption lineshape which is given by the 
Fourier transform of its zero'th moment J{t) = J cr(F,t)(iF. The normalization of cr is thus 
the physically interesting quantity whose Fourier transform gives the spectral lineshape. 

Eq. (jlHj) and Eq. ()18|1 describe very different physical phenomena: The former corresponds 
to an externally driven nonequilibrium system whereas the latter represents a non driven 
system at equilibrium. However, the two equations have a close formal connection. To 
better see the analogy we recast the solution of Eq. (fTH|) as 

(T(F,t) = {exp [-i [ rfrf/(r)]) (19) 
Jo 

where the averaging is over a subensemble of trajectories which assume the same value F at 



time t. For comparison we write the solution of Eq. ()13p as 



—7 / drde/d 



r 



). (20) 



The average here is over a subensemble of trajectories (Fig. 1) which start at thermal 
equilibrium at t = and end up in state m at time t. We reiterate that de/dr only 
corresponds to the continuous part of the trajectory; the jumps should be excluded. 

In both equations we compute a phase space distribution whose time dependent normal- 
ization gives the generating function for the quantity of interest (either work or absorption 
spectrum). The first term in the r.h.s. of Eq. (|13|)) or (|18|) represents a regular dynam- 
ics which conserves the normalization, whereas the second term is responsible for varying 
the normalization, which is the physically interesting property. In Eq. p3|l or Eq. this 
comes from the flow of adiabatic state populations in and out of the system, originating from 
the time dependence of adiabatic energies. In the dephasing problem (Eq. (fTHj) or (HD)), the 
coherence a acquires a random phase (U) which is different for different trajectories and it 
thus constitutes path function. This translates to loss of magnitude (dephasing and deco- 
herence) upon ensemble averaging. The work in the nonadiabatic dynamics is related to an 
integral over the continuous part of the energy trajectory in the same way that the phase 
in the lineshape problem is obtained by an integral over the fluctuating two level frequency. 
This analogy suggests that it should be possible to carry out nonequilibrium lineshape mea- 
surements by e.g. switching on an electric field for a solute in a polar solvent and observing 
the spectroscopic analogues of the Jarzynski relations by looking, for example, at the time 
dependent Stokes Shift. 
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FIG. 1: A "classical trajectory" in population space of adiabatic states. Shown is a trajectory 

which starts at state |(/9„(0)) and ends in state \^Pm{t)) and undergoes j = 5 jumps at time ti t^. 

The solution of the master equation is given by a path integral over all possible trajectories (sum 
over all possible values of j = 0, 1, 2 and integration over tj) 
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